Disordering effects of color in nonequilibrium phase transitions 
induced by multiplicative noise 
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The model introduced by Van den Broeck, Parrondo and Toral [Phys. Rev. Lett. 73, 3395 (1994)] — leading to 
a second-order-like noise-induced nonequilibrium phase transition which shows reentrance as a function of the 
(multiplicative) noise intensity a — is investigated beyond the white-noise assumption. Through a Markovian 
approximation and within a mean-field treatment it is found that — in striking contrast with the usual behavior 
for equilibrium phase transitions — for noise self-correlation time r > 0, the stable phase for (diffusive) spatial 
coupling D — > oo is always the disordered one. Another surprising result is that a large noise "memory" also 
tends to destroy order. These results are supported by numerical simulations. 
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Too often do we resort, in studying nonequilibrium sys- 
tems, to the paradigmatic body we have inherited from equi- 
librium thermodynamics. Though most times this way of rea- 
soning is of valuable help for us to interpret the results, we 
should be more aware of the fact that sometimes it can be seri- 
ously misleading. An archetypical example is the intuitive im- 
age we have developed of a close relationship between noise 
and disorder, and between spatial coupling (and also between 
time-correlation) and order. Regarding the first, whereas it 
is true that studies on e.g. Ginzburg-Landau models subject 
to additive noise seem to reinforce this "rule", [[l|-|3|] in the 
last decade we have also witnessed examples of exactly the 
opposite trend — namely, dynamical systems in which a multi- 
plicative noise couples to the system's nonlinearities in such a 
way that it generates a transition towards an ordered state. In 
fact, it is by now a well-known fact that the noise can induce a 
unimodal-bimodal transition in some O-dimensional models; 
[Q] nevertheless, this result can still be argued to be some- 
what restricted since in this case there cannot be breakdown 
of ergodicity, which is required for a true (nonequilibrium, 
noise-induced) phase transition. 

Recently, a model was introduced whereby an extended sys- 
tem subject to a Gaussian multiplicative noise — white both 
in space and time — can undergo a noise-induced symmetry- 
breaking transition towards an ordered state: this became the 
first example of a purely noise-induced, nonequilibrium, or- 
dering phase transition. [|]| This result was obtained within a 
mean-field approximation and confirmed afterwards through 
extensive simulations in d = 2. [^] In this case, and at vari- 
ance with the case of order-disorder transitions at equilib- 
rium (induced as we know by the spatial coupling constant 
D and the bistability of the local potential) it is the short-time 
instability induced by the multiplicative noise intensity a — 
reinforced by the spatial coupling D — which induces the tran- 
sition. Neither the d = system (D = 0) nor the determin- 
istic one (a = 0) show any transition; moreover — and strik- 
ingly enough — those systems exhibiting noise-induced transi- 



tions in d = are automatically ruled out as candidates for 
this phenomenon. Such a noise-induced phase transition — 
besides being of a second-order type as a function of the noise 
intensity — has the noteworthy feature of being reentrant: the 
ordered state can be found only inside a window determined 
by two values of a. A similar reentrant effect has been ob- 
served in the Ginzburg-Landau model with multiplicative and 
additive noises. 

One can question whether it is realistic enough to con- 
sider a genuine multiplicative noise as white. It appears more 
likely that the kind of fluctuations leading to multiplicative 
noise — coming in general from the coupling with an exter- 
nal source — will exhibit some degree of spatiotemporal cor- 
relation. [p|,p[-^0|] Moreover, one should expect new nontriv- 
ial effects as a consequence of the color of the multiplicative 
noise: after all — whereas it was shown in Refs. [||3|] that an 
ordering nonequilibrium phase transition can be induced in a 
Ginzburg-Landau model by varying the correlation time of the 
additive noise — a reentrant behavior has been found recently 
in a (d = 0) colored-noise-induced transition. JTT| ] 

It is our aim in this work to investigate the effects of 
the self-correlation time r of the multiplicative noise on the 
model of Refs. [pra]. To that end we use an Ornstein- 
Uhlenbeck (OU) noise and apply — in the framework of a 
mean-field treatment — a "unified colored-noise approxima- 
tion" (UCNA), ^ together with an interpolation scheme 
that extends its range of validity in r. [ pj| ] Our main find- 
ing is that — at variance with the usual behavior in equilibrium 
statistical mechanics — a large coupling constant D leads in- 
variably for r > to a disordered state. Since (as discussed 
thoroughly in Ref. |Q|) the clue for the phase transition seems 
to reside in an instability occurring in the short-time behavior, 
and the model introduced in Ref. [|J was precisely chosen as a 
representative (perhaps the simplest one) of a host of systems 
exhibiting such an instability — and hopely, a noise-induced 
phase transition with similar characteristics — this unexpected 
result warns experimentalists seeking for concrete realizations 
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of this phenomenon not to tune naively the spatial coupling in- 
tensity D up to a very large value (as one would do guided by 
"experience") but to look instead for an optimal value of D 
for which the order parameter would take its maximum value. 

Another striking result is that — consistently with the re- 
sult in Ref. ijlljj — increasing the "memory", i.e. the self- 
correlation time t of the noise, does not favor (as one would 
naively expect) the transition towards an ordered phase but 
all the way around. This is indicated by the following facts: 
(a) the threshold (critical) value of a is a strongly increasing 
function of r, and (b) the window in D available to the or- 
dered phase strongly shrinks as r increases. There is nonethe- 
less a hint that (like in Ref. [|TT|]) a small amount of color 
could induce order slightly beyond the upper critical value of 
a corresponding to r = 0, but since the region in which this 
phenomenon occurs is somewhat narrow, and the comparison 
with simulations in d = 2 — together with finite-size scaling — 
made in Refs. [B^] sheds much doubt on the precise location 
of this value, we prefer not to take this result too seriously for 
the moment (although it certainly deserves further analysis). 

As in Ref. [||] we shall resort to a lattice version of the ex- 
tended system, whose state at time t will then be given by 
the set of stochastic variables {xi(t)} (i = 1, . . . , L d ) defined 
at the sites of a hypercubic <i-dimensional lattice of side 
L. The variables {xi} obey the following system of ordinary 
stochastic differential equations (SDE): 

%i = f{xi) + g[xi)r}i + — ( x i ~ Xi ) W 

jen(i) 

where D is the lattice version of the diffusion coefficient, n(i) 
stands for the set of 2d sites which form the immediate neigh- 
borhood of site Yi, and rji is the colored multiplicative noise 
acting on site r^. This coupled set of Langevin-like equa- 
tions is the discrete version of the partial SDE which in the 
continuum would determine the state of the extended system, 
the last term being replaced — in the continuum limit — by the 
Laplacian operator V 2 x. The specific case analyzed in Ref. 
[||] (which the authors conjecture that could be the simplest 
example exhibiting such a transition) is 

f(x) = -x(l + x 2 ) 2 and g(x) = l + x 2 (2) 

As in Refs. [g,|Tl|], the noises {iji} are taken to be OU ones, 
i.e. Gaussian-distributed stochastic variables with zero mean 
and the following correlations: 

(■mWnjV)) = s i ~e W (-^-^) O) 

In the limit t — > the OU noise r\i (i) tends to the white noise 
if{t) with correlations (^(iJfTCO) = o 2 b t] b(t - f), 
which is the case studied in Ref. [pp. 

The non-Markovian character of the process {x{\ (due to 
the colored noise {rji}) makes it difficult to study. However, 
there are some approximate Markovian techniques that — 
whereas capturing some of the essential features of the com- 



plete non-Markovian process — strongly simplify the treat- 
ment of the equations, allowing to exploit well-known Marko- 
vian techniques. jl0| ] Amongst those approximations, the 
UCNA and related interpolation schemes are very useful since 
they can reproduce the limits of small and large correlation 
time r. [0113]] As discussed in Ref. [Q for a single SDE, 
the conditions assumed in the UCNA indicate that its validity 
should decrease with increasing noise intensity. On the other 
hand, regarding the r-dependence, the UCNA becomes exact 
for r — > and for r — > oo. Although the interpolation pro- 
cedure in Ref. jl3| ] extends the validity range of this effective 
Markovian approximation, it is still not clear how far it does 
so. 

We now sketch the main lines of our calculation (a more 
detailed account is given in Ref. 0]): 

i) For our particular problem, the UCNA proceeds by tak- 
ing the time derivative of Eqs.(|l]) and — after substitution of 
the Langevin equations satisfied by the OU noise (r?)j = 
—rji + a^i, where are standard white-noise variables) — 
setting to zero not only '±i (a usual adiabatic elimination) but 
also (ii) 2 , m order to recover a proper Fokker-Planck equa- 
tion (FPE) description. [ |l5| ] 

ii) In order to reduce the complexity of the resulting system 
of Markovian SDE, we make the approximation of replac- 
ing (under the hypothesis that the system is isotropic) in each 
equation of that set the 2c? variables Xj by a single one yi. 

iii) Within the Stratonovich prescription we are left with the 
FPE for a bivariate steady-state probability distribution func- 
tion (pdf) P st (xi,yi). To the drift and diffusion coefficients 
of this FPE we apply an approximation in the spirit of the 
Curie- Weiss mean-field type of approach used in Ref. [B], 
so deriving an effective stationary joint pdf P st (x,y) (we 
have dropped the subindex i for brevity), from which we 
derive a one-site pdf P st (x; (x)) by assuming P st (x,y) = 
P st (x)S(y-(x)). 

iv) The value of (x) follows then from a self-consistency rela- 
tion similar to that of Ref. [^] : 

(4, 

This equation has always the trivial solution (x) = corre- 
sponding to a disordered phase. When other stable, nontrivial, 
(x) ^ solutions appear, the system develops order through a 
genuine phase transition and m = \ (x) \ can be considered as 
the order parameter (due to the symmetry of the problem, both 
±(a;) are solutions of the previous equation). In the white- 
noise limit t = this is known to be the case for sufficiently 
large values of the coupling D and for a window of values of 
the noise intensity a E [<7i, 02]- 

We now discuss how the presence of ordered states is al- 
tered by nonzero values of t in the mean-field study. Figure [l] 
shows, in the parameter subspace a — D, the boundaries sep- 
arating the ordered and disordered phases for different values 
of t. The noteworthy aspects of this graph are the following: 
i) For fixed a > 1 and r > the ordered states can exist only 
within a window of values for D. In other words, the noise- 
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induced nonequilibrium phase transition exhibits reentrance 
not only with respect to a (as in the r = case) but also with 
respect to D. 

ii) For fixed D and a inside the r = phase boundary — as 
indicated, for example, by point (a) in Fig.|l] — there always 
exists a value of the correlation time r beyond which the sys- 
tem becomes disordered. Furthermore, there seems to exist a 
value of t > 0.123 beyond which order is impossible, what- 
ever the values of a and D. 

iii) For fixed (and large enough) values of D, and for values 
of a that would correspond to the disordered phase for r = 0, 
a small increase in r induces a transition towards an ordered 
phase — as indicated by the point marked (b) in Fig.[j]. How- 
ever, a further increase in t can again lead to disorder. In other 
words, the transition can also be reentrant with respect to r. 
Regarding the reentrant nature of the transition with respect 
to D, in Fig.|l]it can be seen that — as r increases from zero — 
the maximum value of D compatible with the ordered phase 
reaches, for a large enough, a "plateau" which is a decreasing 
function of r. At the same time, the minimum value of D (that 
at t = goes like D oc erf) tends also to become constant as 
a function of a as r increases, so shrinking the window avail- 
able for the ordered phase until it virtually disappears. 



course, consistent with what occurs around point (b) in Figjj]. 
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FIG. 2. Mean-field prediction for the order parameter m as a 
function of the spatial coupling D, for noise intensity a = 2 and 
self-correlation times r = (dashed line) and r = 0.01 (continuous 
line). Notice that whereas for r = the curve tends to the asymptotic 
value ( cr 2 — 1 ) 1 / 2 = 1 . 73 , f or t = . 1 the order parameter seems to 
fall off rather abruptly to zero for a large value Dmf of the coupling. 
Simulation results for different system sizes: L = 16 (asterisks), 
L = 32 (rhombs) and L — 64 (triangles) are also included, showing 
that m indeed decays to zero, although much slowly and for a much 
larger "critical" value D c . 
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FIG. 1. Projection of the mean-field phase diagram onto the 
a — D plane, for r = (continous line), r = 0.015 (dotted line), 
r = 0.05 (dashed line), r = 0.1 (dotted-dashed line), and r = 0.123 
(triple dotted-dashed line). For each curve, the ordered zone is the 
area inside the curve (for r = we have marked the ordered and 
disordered regions with "o" and "d", respectively). Points (a) and (b) 
correspond to the transitions referred to in the text. 

Regarding the character of the transition at r = 0, we have 
checked that as r decreases the lower-left elbow climbs up 
the a\ branch corresponding to r = 0, but the slope of the 
corresponding t^O branch is always positive. 

The previous features can also be inferred from the behav- 
ior of the order parameter m. Plotting mvs. cr for a fixed value 
of D and different values of r, we would see a general trend 
of the ordered zone to shrink and disappear with increasing r. 
Whereas the lower critical value <j\ increases monotonically 
with r, the upper value cr% first increases a little and then be- 
comes a monotonically decreasing function of r. This is, of 



Since the previous results have been obtained in the mean- 
field and UCNA approximations, and their range of validity 
is somewhat unclear, we have also performed numerical sim- 
ulations in order to have an independent check of the predic- 
tions. As a representative example (corresponding to the phe- 
nomenon (i) above — namely, the destruction of the ordered 
phase by an increasing coupling constant D) we plot in Fig.|| 
771 vs. D as predicted by our mean-field theory, and results 
coming from a numerical integration of the SDE, for a fixed 
and two values of r. Although for t^O the numerical results 
do not follow the mean-field theory, it is obvious that there is 
an optimal value of the coupling D for which the order pa- 
rameter takes a maximum value, and that order disappears for 
D large enough. From Fig.|| one cannot decide whether the 
maxima will accompany the r = curve as r — > 0. It could 
well be that the phase transition at r = for a fixed and D 
large enough be even a. first-order one. This certainly calls for 
further investigation. 

We stress again the fact that these effects of a colored mul- 
tiplicative noise on an extended dynamical system (unable to 
undergo any phase transition in the absence of noise) are qual- 
itatively different to the ones observed in (nonequilibrium) 
phase transitions driven by a colored additive noise on a proto- 
typic model for equilibrium phase transitions. Whereas 
in the last case the role of the correlation time is to stabi- 
lize the ordered phase and/or induce order in systems that are 
disordered for r = 0, the main effect of color in our case 
is to destroy order. Also, we should not be turned back by 
the quantitative disagreement between the mean-field theory 
and the numerical simulations: it is known that in equilibrium 
phase transitions, mean-field theory overestimates the ordered 
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region and — for example — in the previous study of the same 
model with white noise, [|[|]| the mean-field prediction for the 
upper critical value 02 for the reentrant transition was thrice 
the one found numerically. Although the numerical results 
are affected by finite-size effects — as one would expect in a 
second-order phase transition — one can see unambiguously 
in Fig.|| the decrease of the order parameter with increasing 
coupling D, for t as small as 0.01. 

In order to understand this sudden change in behavior as 
soon as a tiny self-correlation is present, we have studied the 
time evolution equation for (x) (which is small when the pa- 
rameters are around the phase boundary) within the mean- 
field approximation, as ^ — ► 00. In Ref. [|| this simple 
linear (i.e., up to first-order in (x)) criterion of stabilisation 
of the disordered phase was introduced as a way of determin- 
ing the region of appearance of ordered phases. For t <C 1 
and Sz — > 00 (we indeed assume tD and a 2 to be 0(1)), it 
reads 

(x) = -a{x), with a = 1 - ^— °" . (5) 

1 + tD 

When 1 + tD > a 2 it is a > 0, and hence the disordered 
phase ((x) — 0) is stable. On the other hand, if 1 + tD < a 2 
it is a < 0, and it is the ordered phase ( (x) 7^ 0) which be- 
comes stable. In summary, whereas the noise intensity a has 
a stabilizing effect on the ordered phase, as soon as r 7^ the 
spatial coupling D tends to destabilize it. For r = the last 
effect is not present, being then the condition for ordering that 
a > 1 (this is the effect that was reported in Refs. [|]||). Con- 
sidering that the effect of even a tiny correlation is enhanced 
by D, we can understand the abrupt change shown in Figs.|l] 
and H as soon as t 7^ 0. 

This work has focused on the effects of a self-correlation in 
the multiplicative noise on the reentrant noise-induced phase 
transition reported in Ref. [Q]. It appears that for r 7^ 0, a 
strong enough spatial coupling is capable of destroying the 
order established as a consequence of the multiplicative char- 
acter of the noise. The foregoing result can be understood 
by recalling the fact that the ordered phase arises as a conse- 
quence of the collaboration between the multiplicative charac- 
ter of the noise and the presence of spatial coupling. When no 
self-correlation is present, the disordering effect of D cannot 
be felt. This explains the results in Ref. which have been 
rightly interpreted in terms of a "freezing" of the short-time 
behavior by a strong enough spatial coupling. As r increases, 
the minimum value of D required to destabilize the ordered 
phase becomes lower and lower. In this way, the region in 
parameter space available to the ordered phase shrinks further 
and further until it vanishes. 

The main lesson one can draw from the present results is 



that the conceptual inheritance from equilibrium thermody- 
namics (though often useful) is not always applicable. By fol- 
lowing the equilibrium-thermodynamic lore, one would tend 
to think that as D — > 00 an ordered situation is favored. This 
is certainly true for the Curie-Weiss-type models, since in that 
case the deterministic potential is itself bistable and an in- 
crease of spatial coupling has the effect of rising the potential 
barrier between the stable states. In the case we are dealing 
with, the deterministic potential is monostable and it is the 
combined effects of the multiplicative noise and the spatial 
coupling that induce the transition. 
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